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This paper, devoted to the study of spectral pollution, contains both abstract results and applica- 
tions to some self-adjoint operators with a gap in their essential spectrum occuring in Quantum 
Mechanics. 

First we consider Galerkin basis which respect the decomposition of the ambient Hilbert space 
into a direct sum = PS^ ® (1 — P)f), given by a fixed orthogonal projector P, and we localize 
the polluted spectrum exactly. This is followed by applications to periodic Schrodinger operators 
(pollution is absent in a Wannier-type basis), and to Dirac operator (several natural decompositions 
are considered). 

In the second part, we add the constraint that within the Galerkin basis there is a certain 
relation between vectors in Pf) and vectors in (1 — P)f). Abstract results are proved and applied 
to several practical methods like the famous kinetic balance of relativistic Quantum Mechanics. 

© 2008 by the authors. This paper may be reproduced, in its entirety, for non-commercial purposes. 
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Introduction 

This paper is devoted to the study of spectral pollution. This phenomenon of high interest 
occurs when one approximates the spectrum of a (bounded or unbounded) sclf-adjoint op- 
erator A on an infinite-dimensional Hilbert space Sj, using a sequence of finite-dimensional 
spaces. Consider for instance a sequence {Vn} of subspaces of the domain D{A) of A such 
that Vn C Vn+i and Pv„ 1 strongly (we denote by Py^ the orthogonal projector on T^). 
Define the n x n matrices An Pv„APv„. It is well-known that such a Galerkin method 
may in general lead to spurious eigenvalues, i.e. numbers A € R which are limiting points of 
eigenvalues of An but do not belong to a (A). This phenomenon is known to occur in gaps 
of the essential spectrum of A only. 

Spectral pollution is an important issue which arises in many different practical situ- 
ations. It is encountered when approximating the spectrum of perturbations of periodic 
Schrodinger operators [4] or Strum-Liouville operators [Ml |36l [1] . It is a very well reported 
difficulty in Quantum Chemistry and Physics in particular regarding relativistic computa- 
tions [131 dH HH IMl [m HZl 132 ■ It also appears in elasticity, electromagnetism and hydro- 
dynamics; see, e.g. the references in [5]. Eventually, it has raised as well a huge interest in 
the mathematical community, see, e.g., [H [H Hl^l [IS (Ml US] ■ 

In this article we will study spectral pollution from a rather new perspective. Although 
many works focus on how to determine if an approximate eigenvalue is spurious or not (see, 
e.g., the rather successful second-order projection method [231ll])j we will on the contrary 
concentrate on finding conditions on the sequence {Vn} which ensure that there will not be 
any pollution at all, in a given interval of the real line. 

Our work contains two rather different aspects. On the one hand we will establish some 
theoretical results for abstract self-adjoint operators: we characterize exactly (or partially) 
the polluted spectrum under some specific assumptions on the approximation scheme as 
will be explained below. On the other hand we apply these results to two important cases 
of Quantum Physics: perturbations of periodic Schrodinger operators and Dirac operators. 
For Dirac operators, we will show in particular that some very well-known methods used 
by Chemists or Physicists indeed allow to partially avoid spurious eigenvalues in certain 
situations, or at the contrary that they are theoretically of no effect in other cases. 

Let us now summarize our results with some more details. 

Our approach consists in adding some assumptions on the approximating scheme. We 
start by considering in Section [2] a fixed orthogonal projector P acting on the ambiant 
Hilbert space Sj and we define P-spurious eigenvalues A as limiting points obtained by a 
Galerkin-type procedure, in a basis which respects the decomposition associated with P. 
This means A — lim„^oo A„ with A ^ (t{A) and A„ G a{Pv„APY^), where Vn = Vn ® Vn 
for some Vn C Sj^ := PSj and Vn C S)~ := (1 — P)Sj. We show that, contrarily to the 
general case and depending on P, there might exist an interval in M in which there is never 
any pollution occuring. More precisely, we exactly determine the location of the polluted 
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Imposed splitting 
of Hilbert space 


External potential V 


Spurious spectrum 
in the gap (1, 1) 


none 






upper/k 

( fn'^ 

\ . 


jwer spinors 


V ^0 





V bounded 


(-i,-i + sup(y)] 

U[l + inf(y),l) 


unbounded (ex: Coulomb) 


(-1,1) 


dual decomposition iSSf 

\ea-pipnj ' \ Xn J 
< e < 1 


V ^0 





V bounded 


(-l,-2/e+l + sup(F)] 

U[2/e - l + inf(V^),l) 


unbounded (ex: Coulomb) 


(-1,1) 


free decomposition 


any 






Table 1. Summary of our results from Section [2.3l for the Dirac operator + V, when a splitting is imposed 
on the Hilbert space L2(r3,C*). 



spectrum in Section 12.11 and we use this in Section 12.21 to derive a simple criterion on P, 
allowing to completely avoid the appearence of spurious eigenvalues in a gap of the essential 
spectrum of A. 

Then we apply our general result to several practical situations in Section 12.31 We in 
particular show that the usual decomposition into upper and lower spinors a priori always 
leads to pollution for Dirac operators. We also study another decomposition of the ambient 
Hilbert space which was proposed by Shabaev et al [32] and we prove that the set which is 
free from spectral pollution is larger than the one obtained from the simple decomposition 
into upper and lower spinors. Eventually, we prove that choosing the decomposition given 
by the spectral projectors of the free Dirac operator is completely free of pollution. For the 
convenience of the reader, we have summarized all these results in Table [TJ 

As another application we consider in Section 12.3.11 the case of a periodic Schrodinger 
operator which is perturbed by a potential which vanishes at infinity. We prove again that 
choosing a decomposition associated with the unperturbed (periodic) Hamiltonian allows 
to avoid spectral pollution, as was already demonstrated numerically in [6] using Wannier 
functions. 

In Section [3l we come back to the theory of a general operator A and we study another 
method inspired by the ones used in quantum Physics and Chemistry. Namely, additionaly to 
a splitting as explained before, we add the requirement that there is a specific relation (named 
balance condition) between the vectors of and that of Sj~^ . This amounts to choosing a 
fixed operator L : ^+ ?)~ and taking as approximation spaces Vn = ® LV^^ . We do 
not completely characterize theoretically the possible spurious eigenvalues for this kind of 
methods but we give necessary and sufficient conditions which are enough to fully understand 
the case of the Dirac operator. In Quantum Chemistry and Physics the main method is the 
so-called kinetic balance which consists in choosing L — (t(— iV) and the decomposition into 
upper and lower spinors. We show in Section [3. 2. II that this method allows to avoid spectral 
pollution in the upper part of the spectrum only for bounded potentials and that it does 
not help for unbounded functions like the Coulomb potential. We prove in Section [3.2.2l that 
the so-called (more complicated) atomic balance indeed allows to solve this problem also for 
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Balance 
condition 


External potential V 


Spurious spectrum 
in the gap (1, 1) 


kinetic balance 

( \ 
V / \(^-PVn) 


V bounded with 
-1 + sup(y) < 1 + M{V) 


(-1,-1 + supCl/)] 


< K < V3/2 


(-1,1) 


atomic balance 
(^0 )' (^^a-pv'j 


V such that 
— jfj < Vix) where 
< K < V3/2, 
and sup(y) < 2 


(-l,-l + sup(-l/)] 


dual kinetic balance [32j 
\ea ■ pifn J ' \ ipn J 


V bounded 


(-l,-2/e+l + sup(F)] 

U[2/e - l + inf(l/),l) 


unbounded (ex: Coulomb) 


(-1,1) 



Table 2. Summary of our results for the Dirac operator D'^ + V when a balance is imposed between vectors 
of the basis. 



Coulomb potentials, as was already suspected in the literature. Eventually, we show that the 
dual kinetic balance method of |32| is not better than the one which is obtained by imposing 
a splitting without a priori adding a balance condition. Our results for balanced methods 
for Dirac operators are summarized in Tabled 

We have tried to make our results sufficiently general that they could be applied to other 
situations in which there is a natural way (in the numerical sense) to split the ambiant 
Hilbert space in a direct sum Sj ~ S)~^ (BSj~. We hope that our results will provide some new 
insight on the spectral pollution issue. 

Acknowledgements. The authors would like to thank Lyonell Boulton and Nabile Boussaid for 
interesting discussions and comments. The authors have been supported by the ANR project AC- 
CQuaRel of the french ministry of research. 

1. Spectral pollution 

In this first section, we recall the definition of spectral pollution and give some properties 
which will be used in the rest of the paper. Most of the material of this section is rather 
well-known [M El [23l [9] . 

In the whole paper we consider a self-adjoint operator A acting on a separable Hilbert 
space Sj, with dense domain D{A). 

Notation. For any finite-dimensional subspace V C D{A), we denote by Py the orthogonal 
projector onto V and by A\y the self-adjoint operator V V which is just the restriction 
to V of PyAPv. 

As A is by assumption a self-adjoint operator, it is closed, i.e. the graph G{A) C D{A) x S) 
is closed. This induces a norm |H|£)(-^-) on D{A) for which D{A) is closed. For any K C D{A), 

we will use the notation K'^^^^ to denote the closure of K for the norm associated with the 
graph of A, in D{A). On the other hand we simply denote by K the closure for the norm of 
the ambient space 



Spectral Pollution and How to Avoid It 



5 



We use like in [23j the notation (Toss(^) to denote the essential spectrum of A union 
—oo (and/or +00) if there exists a sequence of cr{A) 9 A„ — > —00 (and/or +cxd). Finally, 
we denote by Conv(X) the convex hull of any set X C M and we use the convention that 
[c, d] = if d < c. 

Definition 1.1 (Spurious eigenvalues). We say that X G M. is a spurious eigenvalue of 
the operator A if there exists a sequence of finite dimensional spaces {Vn}n>i with Vn C D{A) 
and Vn C Vn+i for any n, such that 

(i) u;:ilK''^^^ - D{A); 
(a) lim dist (A, a{A\y )) — 0; 

(Hi) A a{A). 

We denote by Spu(A) the set of spurious eigenvalues of A. 

If needed, we shall say that X is a spurious eigenvalue of A with respect to {Vn} to 
further indicate a sequence {Vn} for which the above properties hold true. Note that (i) in 
Definition 11.11 implies in particular that we have Un>iVn = f) since D{A) is dense in Sj by 
assumption. 

Remark 1.1. As the matrix of A in a finite- dimensional space only involves the quadratic 
form associated with A, it is possible to define spurious eigenvalues by assuming only that Vn 
is contained in the form domain of A. Generalizing our results to quadratic forms formalism 
is certainly technical, although being actually useful in some cases (Finite Element Methods 
are usually expressed in this formalism). We shall only consider the simpler case for which 
Vn C D(A) for convenience. 

Remark 1.2. If X is a spurious eigenvalue of A with respect to {Vn} and if B — A is 
compact, then X is either a spurious eigenvalue of B in {Vn} or A G a^isciB). One may 
think that the same holds when B — A is only A- compact, but this is actually not true , as 
we shall illustrate below in Remark \2. 7[ 

Remark 1.3. In this paper we concentrate our efforts on the spectral pollution issue, and 
we do not study how well the spectrum a{A) of A is approximated by the discretized spectra 
a{A\Y^). Let us only mention that for every X G <^{A), we have dist(A, (T(v4|y^)) — ^ as 

n 00, provided that U„>iV^ ' — D(A) as required in Definition ] 

The following lemma will be very useful in the sequel. 

Lemma 1.1 (Weyl sequences). Assume that X is a spurious eigenvalue of A in {Vn} as 
above. Then there exists a sequence {xn}n>i C D{A) with a;„ G Vn for any n>l, such that 

(1) Pv„{A — X)xn ~* strongly in Sj; 

(2) \xn \ = 1 for all n > 1; 

(3) x„ ^ weakly in f). 

Proof. It is partly contained in [9]. Let A G Spu(A) and consider Xn & Vn \ {0} C D{A) 
such that Pv„iA — A„)a;„ = with lim„^oo -^n = ^- Dividing by ||a:„|| if necessary, we may 
assume that ||a;„|| = 1 for all n in which case Pv„{A — X)xn — > strongly. As {x„} is 
bounded, extracting a subsequence if necessary we may assume that Xn ^ x weakly in ^. 
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What remains to be proven is that x = 0. Let y e L}m>iVm- Taking n large enough we may 
assume that y ^ Vn- Next we compute the following scalar product 

0= lim {Pv„{A- X)xn,y) = lim {x^, {A ~ X)y) ^ {x, (A - X)y) . 

n — *oo n — >oo 

As Um>iKn is dense in D{A) for the norm of G{A), we deduce that {x, {A — X)y) — for 
all y £ D{A). Thus x £ D{A*) — D{A) and it satisfies Ax = Aa;. Hence x = since A is not 
an eigenvalue of A by assumption. □ 

The next lemma will be useful to identify points in Spu(74). 

Lemma 1.2. Assume that A is as above. Let be an orthonormal system of K 

vectors in D{A) such that xl^ Q for all j = I..K. Denote by Wn the space spanned by 
xl^, x^ . If X £ M. is such that lim„^oo dist (A , cr(A|vi/^)) = 0, then X e Spu(A) U o'{A). 



Proof. Consider any nondecreasing sequence {Vn} such that U„>iK ^ ' = D{A). Next 
we introduce V{ := Vi, mi — and we construct by induction a new sequence {V^} and an 
increasing sequence {to„} as follows. Assume that and m„ are defined. As x'^ for all 
k = we have limm^oo {Ay,x'^) — for all y G and all k = 1..K. Hence the matrix 
of A in 1^ + Wm becomes diagonal by blocks as m cx3. Therefore there exists run+i > nrin 
such that the matrix of A in := V^,' + Wm^+i has an eigenvalue which is at a distance 
< 1/n from A. As 14 C V,' for all n, we have U„>iV^^^'^'' — D{A). By construction we also 
have lim„^oo dist (A , a{A\v^)) — 0. Hence either A e <y{A), or A S Spu(yl). □ 

In the following we shall only be interested in the spurious eigenvalues of A lying in the 
convex hull of acssiA). This is justified by the following simple result which tells us that 
pollution cannot occur below or above the essential spectrum. 

Lemma 1.3. Let X be a spurious eigenvalue of the self-adjoint operator A. Then one has 

Tr(x(-oo,A](^)) =Tr(x[A,oo)(^)) =+oo. (1.1) 
Saying differently, X G Conv ((Tcss(^))- 

Proof. Assume for instance P := X(-oo,a](^) is finite-rank. As A ^ '^{A), we must have 
P = X(-oo,A+e] (A) for some e > 0. Let {xn} be as in Lemma fTTT] As P is finite rank, Pxn — + 
and {A — X)Pxn strongly in Sj. Therefore ^V„(^ ~ A)P^x„ strongly. Note that 
(A-A)P-L > eP-L, hence (A - A)P^a;„, x„) = (P^(A - A)P-^x„, x„) > e||P-^a;„f. 
As the left hand side converges to zero, we infer which contradicts Lemma [TTTl □ 

We have seen that pollution can only occur in the convex hull of acss{A). Levitin and 
Shargorodsky have shown in that (|l.ip is indeed necessary and sufficient. 

Theorem 1.1 (Pollution in all spectral gaps [23j). Let A be a self-adjoint operator on 
^ with dense domain D{A). Then 

Spu(A) U ^ CSS 

{A) = Conv ( (A)). 

Remark 1.4. As J := Conv ((Toss(^)) \ <5'css(^) only contains discrete spectrum by as- 
sumption, Theorem \1.1\ says that all points but a countable set in J are potential spurious 
eigenvalues. 
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/ 

Spii(A) 



Fig. 1. For an operator A which has a spectral gap [a, 6] in its essential spectrum, pollution can occur in the 
whole gap. 



Remark 1.5. It is easy to construct a sequence Vn like in Definition \1.1\ such that 
dist {X,cr{AivJ) for all A e Conv ((7css(^)), see fMI- 

Theorem 11.11 was proved for bounded self-adjoint operators in [55] and generalized to 
bounded non self-adjoint operators in 10\. For the convenience of the reader, we give a 
short 

Proof. Let A e Conv ((Tcss(^)) \ <5'oss(^) and fix some a < A and 6 > A such that a,5 €E 
(Tcss(^) priori we might have b = +oo or a ~ — oo). Let us consider two sequences 
{x„}, {un} C D{A) such that (A - a„)x„ — > 0, (A - 6„)2/„ 0, ||a;„|| = ||?/„|| = 1, a;„ ^ 0, 
2/„ ^ 0, a„ — > a and &„ b. Extracting subsequences if necessary we may assume that 
(xnyVn) ^ as n ^ oo. Next we consider the sequence z„{0) := cos9 Xn + smO j/„ which 
satisfies ||2:n(^)|| ^ 1 and Zn{9) uniformly in 6. We note that {Azn{Q)^ 2n(0)) = a„ + o(l) 
and (Az„(7r/2), z„(7r/2)) = 6„ + o(l). Hence for n large enough there exists a 0„ G (0,7r/2) 
such that (Az„(0„), Zn{9n)) = A. The rest follows from Lemma FOl □ 



2. Pollution associated with a splitting of 

As we have recalled in the previous section, the union of the essential spectrum and (the 
closure of) the polluted spectrum is always an interval: it is simply the convex hull of 
dess{A). It was also shown in [23] that it is possible to construct one sequence {Vn\ such 
that all possible points in Spu(A) are indeed {V„}-spurious eigenvalues. But of course, not 
all {Vn] will produce pollution. If for instance Pv„ commutes with A for all n > 1, then 
pollution will not occur as is obviously seen from Lemma 11.11 The purpose of this section 
is to study spectral pollution if we add some assumptions on {T/J^}. More precisely we will 
fix an orthogonal projector P acting on and we will add the natural assumption that Py^ 
commute with P for all n, i.e. that Vn only contains vectors from Pio and (1 — P)fi. 

As we will see, under this new assumption the polluted spectrum (union (Tcss(^)) will in 
general be the union of two intervals. Saying differently, by adding such an assumption on 
{Ki}, we can create a hole in the polluted spectrum. A typical situation is when our operator 
A has a gap in its essential spectrum. Then we will see that it is possible to give very simple 
condition^ on P which allow to completely avoid pollution in the gap. 

Note that our results of this section can easily be generalized to the case of a partition of 
unity {PiYi^i of commuting projectors such that 1 = Y^^i=i Pi- Adding the assumption that 
Pv^ commutes with all P^'s, we would create p holes in the polluted spectrum. This might 
be useful if one wants to avoid spectral pollution in several gaps at the same time. 



^Loosely speaking it must not be too far from the spectral projector associated with the part of the spectrum 
above the gap, as we will see below. 
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2.1. A general result 

We start by defining properly P-spurious eigenvalues. 

Definition 2.1 (Spurious eigenvalues associated w^ith a splitting). Consider an or- 
thogonal projection P : ^ f). We say that A G M is a P-spurious eigenvalue of the 
operator A if there exist two sequences of finite dimensional spaces {V^,^}ri>i C PS) H D{A) 
and {V^}n>i C (1 — P)^ n D{A) with C for any n, such that 

l—DiA) 

(1) yJn>i{Vn ®V+) ^D{A); 

(2) £m dist (\,a(A^^y+^y-^)) = 0; 

(3) \ia{A). 

We denote by Spu(yl, P) the set of P-spurious eigenvalues of the operator A. 

Now we will show as announced that contrarily to Spu(^) U (Tcss(^) which is always an 
interval, Spu(A, P) U (Toss(^) is the union of two intervals, hence it may have a "hole". 

Theorem 2.1 (Characterization of P-spurious eigenvalues). Let A be a self-adjoint 
operator with dense domain D(A). Let P be an orthogonal projector on f) such that PC C 
D[A) for some C C D{A) which is a core for A. We assume that PAP (resp. (1 — P)^(l — 
P)) is essentially self-adjoint on PC (resp. (1 — P)C), with closure denoted as A^p^^ (resp. 
A\(i_P'^^). We assume also that 

inf (Tess(^l(i-P)ft) < inf (7oss(^|p^,)- (2.1) 

Then we have 

Spu(A, P) U (Tcss(A) = [inf (Toss(^),SUp(Toss(^|(l-P)i3)] 

U [inf (Toss (^|P^), sup (Tess(^)] ■ (2.2) 



Spu(A,P) 




Fig. 2. Illustration of Theorem 12.11 for an operator A with a gap [a, b] in its essential spectrum, pollution 
can occur in the whole gap, except between the convex hulls of (Tess(^|pj5) and iTess (^| (i-p)j5 ) . 



Let us emphasize that condition (|2.ip always holds true, exchanging P and 1 — P if 
necessary. Usually we will assume for convenience that 1 — P is "associated with the lowest 
part of the spectrum" in the sense of (|2.ip . 

As mentioned before, an interesting example is when A possesses a gap [a, b] in its 
essential spectrum, i.e. such that (a, b) n f7oss(^) = and 

Tr {x{-oo,a]{A)) = Tr(x[b,oo)(^)) = +oo. 

Then taking 11 = X[c,oa){A) and C = D{A) we easily see that Spu(^,n) n (a, 6) = 0. The 
idea that we shall pursue in the next section is simply that if P is "not too far from 11", 
then we may be able to avoid completely pollution in the gap [a, 5]. 
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Before writing the proof of Theorem 12. li we make some remarks. 

Remark 2.1. If the symmetric operators PAP and (I — P)A{1 — P) are both semi-bounded 
on their respective domains PC and (1 — P)C, then the inclusion C in p.2p is also true 
provided that A\p^ and A^^i^p-^fj are defined as the corresponding Friedrichs extensions. 
The essential self-adjointness is only used to show the converse inclusion D. 

Remark 2.2. An interesting consequence of Theorem ] 2. 1\ is that the set of spurious eigen- 
values varies continuously when the projector P is changed (in an appropriate norm for 
which the spectra of A^p^j and A|(i_p)^ change continuously). This has important practi- 
cal consequences: even if one knows a projector which does not create pollution, it could in 
principle be difficult to numerically build a basis respecting the .splitting of Sj induced by P. 
However we know that pollution will only appear at the edges of the gap if the elements of 
the Galerkin basis are only known approximately. 

Proof. We will make use of the following result, whose proof will be omitted (it is an 
obvious adaptation of the proof of Lemma II. 2p : 

Lemma 2.1. Assume that A is as above. Let (x,;',, ...,x^) and (y^, ) be two orthonor- 

mal system^ in PS)nD{A) and (1 — P)Sjr\D{A) respectively, such that x^^ ^ and y^ ^ 
for all j — 1..K and j' = 1..K' . Denote by Wn the space spanned by x\, ...,x^ ,yl^, ...,y^ . 
If X eR is such that lim„_»oo dist (A , a{A\w„)) = 0, then X e Spu(A, P) U cr(A). 

In the rest of the proof, we denote [a,b] :— Conv {acss{A)) , [ci,di] := 
Conv ((Jcss(^l(i-p)fi)) and [c2,c?2] := Conv ((Tess(^|pfl))- For simplicity we also introduce 
c = min(ci, C2) ~ ci, and d — max(c?i, d2). Recall that we have assumed ci < C2. 

Step 1. First we collect some easy facts. The first is to note that Spu(yl, P) C Spu(A) C 
[a, 6], where we have used Theorem ll.il Next we claim that 

[ci, di] U [c2, da] C Spu(A, P) U a{A) n [a, b]. (2.3) 

This is indeed an obvious consequence of Thcorem ll.il applied to A^pf, and A|(i_p)f|, and 
of Lemma 12.11 

Step 2. The second step is less obvious, it consists in proving that 

[a, c] U [d, 6] C Spu(yl, P) U a{A) n [a, b] (2.4) 

which then clearly implies 

[a, di] U [c2, 6] C Spu(A, P) U a{A) n [a, b]. 

Let us assume for instance that d < b and prove the statement for [d, b] (the proof is 
the same for [a, c]). Note that b may a priori be equal to +00 but of course we always have 
under this assumption d < +00. In principle we could however have d = —00. In the rest of 
the proof of (|2.4p . we fix some finite A e {d,b) and prove that A G Spu(A, P) U cr{A). We 
also fix some finite d' such that d < d' < X. We will use the following 



''We will allow K = or K' = 0. 
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Lemma 2.2. Assume that b G (Tcss(^)- Then there exists a Weyl sequence {xn} C C such 
that {A ~ bn)xn 0, ||x„|| — 1, x„ ^ 0, bn ^ b and 

^ and „i^~^j^"„ - weakly. (2.5) 



||(l-P)a;„|| 

Proof. Let bn b and {yn} C C be a Weyl sequence such that {A — hn)yn ^ with 
||y„|| = 1, yn (note we may assume {y„} C C since C is a core for A). We denote 
Vn = ut + Vn where ?/+ G PC C D{A) and j/^ G PC C D{A). Extracting a subsequence, 
we may assume that ||y^||^ and that ||y^||^ — > t~\ note ^+ + = 1. It is clear that if 

l'^ > 0, then y^ \\yt\\~^ ~^ since 0. We will assume for instance f + = and £^ = 1. 

Next we fix an orthonormal basis {e^} C PC of P9j, we define 

k 

'■— {ei,ynk)Si 
i=l 

and note that 

fc 

(A - 6„Jr+ = ^ (^(e„y„,)Ae, + (e„ (A - 6„ J?/„ Je, - (Aej,y„Je,y 
1=1 

For /c fixed and any i = l..fc, we have 

lim (ei,y„)= lim (e^, (A - 6„)2/„) lim (Aei,2/„)=0. 

n — *oo n — >oo n — >oc 

Hence, for a correctly chosen subsequence {y^j.}, we may assume that 
satisfies lim Ik^ll = lim ||(^ — &nfc)?'fc^|| =0. 

Next we define Xk ■= yn^ ~ — iVuk ^ ^t) ^ Vuk which satisfies = 1 + o(l) since 
Ikfe"! ~* ^- construction, we have — — r^l" G span(ei, ...,6^)-'-, hence necessarily 
^k Ikfc II ^ ^ 0. Eventually, we have [A — bn^)xk — > strongly, by construction of r^. □ 

In the rest of the proof we choose a sequence {a;„} like in Lemma 12.21 and denote x^ = 
Pxn and x^ = [1 — P)xn- By the definition of d and the fact that v4|(i_p)f, is essentially 
selfadjoint on (1 — P)C, we can choose a Weyl sequence {y^} C (1 — P)C such that (1 — 
P){A — dn)yn ~^ 0, ||y~|| = 1, weakly and dn di < d. Extracting a subsequence 

from {y~} we may also assume that y~ satisfies 

lim l^^.Vn) = lim IP^^^Vn) = lim /^,y-\ =0 (2.6) 



"Let us now introduce the following orthonormal system 



sm« y„ 



^,z;„(0)l with v^{e):=^= (2.7) 

J 1 + 23fJ cos 6* sin e*/ ||%|| , y„ 



and denote by An{9) the 2x2 matrix of A in this basis, with eigenvalues A„(6') < /i„(6'). As 



INn II ^ ^ weakly, we have 



limsup sup Xn{e) < lim sup ^^^"''^^"^ < d2 < d. (2.8) 

n— >oo ee[0,7r/2] n^oo ||xn|| 



Spectral Pollution and How to Avoid It 



11 



When 6* = 0, wc know by construction of Xn that A„(0) has an eigenvalue which converges 
to & as n — > oo. Since 6 > d by assumption, this shows by (j2.8p that this eigenvalue must 
be /^ri(0), hence we have /l(„(0) ^ 6 as n — > oo. On the other hand, the largest eigenvalue of 
A„(7r/2) satisfies for n large enough 



M„(V2) < max ( ^4r^^ , {Ay-,y-)] + 



W^n I 




< d', 



where we have used (|2.6|) . ||x+|| 0, ^ 0, and the definition of d' > d. 

By continuity of ^in{0), there exists a 0„ e (0,7r/2) such that fj,n{0n) — A. Next we note 
that the two elements of the basis defined in (|2.7p both go weakly to zero by the construction 
of and of y~ . Hence our statement A £ Spu(^, P) U a{A) follows from Lemma \TJ\ 

Step 3. The last step is to prove that when di < C2, 

(di, C2) n (Spu(A, P) U aess(A)) = 

(there is nothing else to prove when C2 < di). We will prove that (rfi, C2) n Spu(74, P) = 0, 
the proof for cross(^) being similar. Note that under our assumption di < C2, we must have 
rfi < 00 and C2 > —00, hence A^p^ and v4|(]^_p-)^ are semi-bounded operators. As noticed in 
Remark l2.1[ it is sufficient to assume for this step that A^pgj and A|(i_p-)^ are the Friedrichs 
extensions of {PAP, PC) and ((1 — P)A{1 — P), (1 — P)C) without assuming a priori that 
they are essentially self-adjoint. 

Now we argue by contradiction and assume that there exists a Weyl sequence S 
C D{A) like in Lemma [TTTl for some A G (rfi, C2). We will write Xn = a;+ with 
^ s-iid ^ • Wc have P^y+^y- (A — X)xn 0, hence taking the scalar product 
with xt and x~ , we obtain 



lim {{A- X)xn,x+) = lim {{A ^ X)xn, x') = 0. (2.9) 

n — >oo n — >oo 

The space C being a core for A, it is clear that we may assume further that x„ G C and still 
that p.9p holds true. In this case we have x^,x^ E D{A) hence we are allowed to write 

((A - A)x+,x+) + {{A - X)x-,x+) ^ 0, 



{{A - A)a;-, a;„ ) + {{A - A).t„ ,x+)^0. 

Taking the complex conjugate of the second line (the first term is real since A is self-adjoint) 
and subtracting the two quantities, we infer that 

{{A - X)x+,x+) - {{A - X)x-,x-) ^ 0. (2.10) 

As by assumption A £ {di, C2), we have as quadratic forms on PC and (1 — P)C, P{A~X)P > 
eP — r and — (1 — P)(A — A)(l — P) > e(l — P) — r' for some finite-rank operators r and r' 
and some e > small enough. Hence we have 

{{A~X)x+,x+) - {{A - X)x-,x-) > e f + e fa:;- f + o(l). 

This shows that we must have a;„ — > which is a contradiction. □ 
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2.2. A simple criterion of no pollution 

Here we give a very intuitive condition allowing to avoid pollution in a gap. 

Theorem 2.2 (Compact perturbations of spectral projector do not pollute). Let 

A be a self-adjoint operator defined on a dense domain D{A), and let a < b be such that 

(a,6)nacss(A) =0 and Tr (x(_oo,a] (^)) = Tr (x[b,oo)(^)) = +oo. (2.11) 

Let c G (a, b) \ cr(A) and denote 11 :— X(c,oo)(^)- Let P be an orthogonal projector satisfying 
the assumptions of Theorem \2.1l We furthermore assume that {P -n)\A- c\^/^, initially 
defined on D{\A — c|^/^), extends to a compact operator on Sj. Then we have 

Spu{A, P) n (a, b) = 0. 

As we will see in Corollarv l2.1[ Theorem 12.21 is useful when our operator takes the form 
A + B where B is A-compact. Using the spectral projector P = 11 of A will then avoid 
pollution for A + B, when A is bounded from below. 

Remark 2.3. We give an example showing that the power 1/2 in \A — c|^/^ is sharp. Con- 
sider for instance an orthonormal basis {e^} of a separable Hilbert space S), and define 
^ ■= E„>i"|e^)(e,tl- Choosing c = 1/2, we get U ^ X[i/2,oo)(-4) = E„>i |e^>(e,tl- Define 
now a new basis by /+ — cos0„e+ + sin^?„e^, f^ = sin0„e+ — cos9„e~, and introduce 
the associated projector P = J2n>i \fn){fn\- Consider then Vn ■= span{f^, • • • , /^} 
for which we have a{A^Vn) = {0,1,--- ,n — l,nsin^0„}. On the other hand it is easily 
checked that (P — n)|yl — 1/2|" is compact if and only if n°'9n — s- as n —> oo. Hence, if 
< a < 1/2 we can take 9n = l/\/2n and we will have a polluted eigenvalue at 1/2 whereas 
{P - n)|A - 1/2|" is compact. 

We now write the proof of Theorem 12.21 



Proof. We will prove that iTcss(^|pfi) C [6, oo). This will end the proof, by Theorem 12.11 
and a similar argument for A|(i_p)^. Assume on the contrary that A G (— oo, b)C\acss{A\pf,). 
Without any loss of generality, we may assume that c > A (changing c if necessary) . As PC is 
a core for A|p^, there exists a sequence {a;„} C PC such that a;„ ^ weakly in ij, ||a;„|| = 1 
and P{A — A)a;„ strongly in f). We have 



((P - n)(A - A)(P - n)x„, x„) + 25R((P - n)(A - c)na;„, x„) 

+ (n(A - A)nx„, a;„) = {P{A - A)a;„, x„) + (A - c)25R(n.T„, (P - n)a;„ 



(2.12) 



where we note that Pxn — Xn ^ D{A) and Ilxn G D{A) since 11 stabilizes D{A). As 
c ^ o'(A), we have that \A— c\~^^^ is bounded, hence P— 11 must be a compact operator, i.e. 
the last term of the right hand side of (|2.12[) tends to as n ^ oo. By the Cauchy-Schwarz 
inequality we have 



\{{P-U){A-c)Ux„ 



< 



A-c\^/^Ilx„ |A-c|i/2(P-n)a;„ 



(2.13) 



As by assumption (P — n)|yl — c|^/^ is compact, we have that (P — Il){A — A)(P — H) and 
\A — c|^/^(P — n) are also compact operators. Hence 



lim 

n—*oo 



\A 



11/2^ 



(P-H)a;„ = lim ((P - H)(A - A)(P - H)x„, x„) = 0. 
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On the other hand we have U{A - A)n = U{A - c)n + (c - A)n > I1\A - c\U since we have 
chosen c in such a way that c > A, and by the definition of H. Hence by p.l2p we have an 
inequahty of the form 



\A 



2e„ 



\A-c\'^^nx„ 



where hm„_^oo Cn 



hm„ 



0. This clearly shows that 



lim 



\A-c\^^^IlXr^ 



= 0. 



Therefore we deduce nx„ — > strongly, \A — c|^/^ being invertible. Hence 
(P — Tl)xn + Tlxn — > and we have reached a contradiction. 



PXr^ 



□ 



We now give a simple application of the above result. 

Corollary 2.1. Let A be a bounded-below self-adjoint operator defined on a dense domain 
D{A), and let a < b be such that 

(a,5)nt7ess(A) =0 and Tr (x(-oo,a] (^)) = Tr (x[fc,oo)(^)) = (2.14) 

Let c € (a, b) be such that c ^ o'{A) and denote H :~ X(c,oo)(^)- 

Let B be a symmetric operator such that A + B is self- adjoint on D{A) and such that 
(^{A-\- B ^ i)^^ ~ {A — i)~^^\A — c]^^"^ , initially defined on D{\A — c\^^'^), extends to a compact 
operator on Sj. Then we have 

Spu(A + B,n) n {a,b) = 0. 

Proof. Under our assumption we have that {A -\- B — i)~^ — (A — i)~^ is compact, hence 
(Tess + B) = (Tess (A) by Wcyl's Theorem [301 ISj and A -\- B is also bounded from below. 
Changing c if necessary we may assume that c ^ a{A -\- B)U a{A). Next we take a curve 
in the complex plane enclosing the whole spectrum of A and A + B below c (i.e. intersecting 
the real axis only at c and c' < inf (t(A) U (j{A + B)). In this case, we have by Cauchy's 
formula and the resolvent identity 



n-x[c,oo)(A + B))|A- 



|l/2 



1 

2i7r 



1 

A + B 



1 



A + B 



A + B 



A + B 



A-z 
1 



A 



\A- 



\A 



c\"^dz 
A 



|l/2. 



A 



-dz 



Since ^ is bounded (we use here that A is bounded-below), we easily deduce that the above 
operator is compact, hence the result follows from Theorem 12.21 □ 



Remark 2.4. Again the power 1/2 in \A — c|^/^ is optimal, as seen by taking B = —A + 
J2n ''^Ifn) {fn\ whcrc A, /+ and On are chosen as in Remark \2.S\ and V-a '■— {ef , e^_]^, e^}. 

Remark 2.5. Corollarv \2.1\ is a priori wrong when A is not semi-bounded. This is seen by 
taking for instance A ^ E„>i "|e^>(ert|-E„>i "|e;:)(e;:l and B = -A+X;„>i "1/^) 
En>i"l/n )(/„"! where f+ = e+/V2 + e-/V2 and f- = -e-/V2 + e+/V2. A short 
calculation shows that (^{A + B)^^ — A^^^\A\°^ is compact for all Q < a < 1 whereas 
G Spu(yl + i3,n) which is seen by choosing again Vn — {ej*^, e^_]^, e^}. 
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2.3. Applications 

2.3.1. Periodic Schrddinger operators in Wannier basis 

In this section, we show that approximating eigenvalues in gaps of periodic Schrodinger 
operators using a so-called Wannier basis does not yield any spectral pollution. This method 
was already successfully applied in dimension 1 in 6J for a nonlinear model introduced in 
[6] . For references on pollution in this setting, we refer for example to [4] . 
Consider d linearly independent vectors ai, in R'^ and denote by 

if aiZe • • • e fldZ 

the associated lattice. We also define the dual lattice 

^* a*Z e • • • ® a*{L with (a.^, a*) = (27r)%. 

Finally, the Brillouin zone is defined by 

.X e M'' I Ibll = inf lb - fcll 

Next we fix an .if-periodic potential Vpcr, i.e. V^c^ix + a) = V^c^x) for all a E Jff. We will 
assume as usual [30] that 

r p = 2 if rf < 3, 
Vpor e LP{J^) where lp>2 if = 4, 
[p = d/2 if d>5. 

In this case it is known [3D] that the operator 

Ap,, - -A + Vp,r (2.15) 
is self-adjoint on i/^(M'^). One has the Bloch-Floquet decomposition 

where Apor(0 is for almost all ^ G a self-adjoint operator acting on the space 

L| = {w e iioc(K^) I u{x + a) = e-"-«u(a;), Va G ^} . 

For any ^, the spectrum of Apor(^) is composed of a (nondecreasing) sequence of eigenvalues 
of finite multiplicity Xk{0 /* hence the spectrum 



Cr(Apcr) = cross (Apcr) = |J Afc (^) 



fc>l 

is composed of bands. The eigenvalues Afe(^) are known to be real-analytic in any fixed 
direction when Vper is smooth enough [39|, 130] , in which case the spectrum of Apa- is purely 
absolutely continuous. 

The operator (|2.15p may be used to describe quantum electrons in a crystal. It appears 
naturally for noninteracting systems in which case Vpcr is the periodic Coulomb potential 
induced by the nuclei of the crystal. However operators of the form p.lSp also appear in 
nonlinear models taking into account the interaction between the electrons. In this case, the 
potential Vper contains an additional effective (mean-field) potential induced by the electrons 
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themselves [7, 6 . In the presence of an impurity in the crystal, one is led to consider an 
operator of the form 

A^-A + Vper + W. (2.16) 

We will assume in the following that 

W eLP{R'^) + L^{R'^) for some p > max(d/3, 1) 

in which case (Ap^r + W — — {Apcr — i)~^ is (1 — A)^^/^-compact as seen by the resolvent 
expansion [30], and one has 

Cross(^) = a{Apcr). 

However eigenvalues may appear between the bands. Intuitively, they correspond to bound 
states of electrons (or holes) in presence of the defect. By Theorem 1 1.1[ their computation 
may lead to pollution. For a finite elements- type basis, spectral pollution was studied in [4] . 

Using the Bloch-Floquet decomposition, a spectral decomposition of the reference pe- 
riodic operator Ap^r is easily accessible numerically. This decomposition can be used as a 
starting point to avoid pollution for the perturbed operator A. For simplicity we shall as- 
sume that the spectral decomposition of Ap^r is known exactly. More precisely we make the 
assumption that there is a gap between the fcth and the (fc -I- l)st band: 

a supAfc(^) < inf Afc+i(^) := b 

and that the associated spectral projector 

-fpcr ■ — X( — oo.c) (^pcr)j C 

is known. The interest of this approach is the following 

Theorem 2.3 (No pollution for periodic Schrodinger operators). We assume Vper 
and W are as before. Then we have 

Spu(A fpcr) n (a, 6) = 0. (2.17) 

Proof. This is a simple application of Corollarv l2.1l □ 

It was noticed in [6] that a very natural basis respecting the decomposition associated 
with Ppei- is given by a so-called Wannier basis [10]. Wannier functions {wk} are defined 
in such a way that Wk belongs to the spectral subspace associated with the fcth band and 
{wk{- — a)}ae^ forms a basis of this spectral subspace. One can take 

^fc(x)-^£7/fc(^,x)de (2.18) 

where Uk{^,-) G is for any ^ G ^ an eigenvector of Aper(C) corresponding to the fcth 
eigenvalue Afc(^). The so-defined {wk{- — a)\a^^ are mutually orthogonal. Formula (|2.18p 
does not define Wk uniquely since the Mfe(^, a;) are in the best case only known up to a phase. 
Choosing the right phase, one can prove that when the fcth band is isolated from other 
bands, Wk decays exponentially [25] . 

More generally, instead of using only one band (i.e. one eigenfunction Ufc(^, x)), one can 
use K different bands for which it is possible to construct K exponentially localized Wannier 



a^b 
2 
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functions as soon as the union of the K bands is isolated from the rest of the spectrum [26l 
[5]. The union of the K bands is caUed a composite band. 

In our case we typically have a natural composite band corresponding to the spectrum 
of Apci which is below c, and another one corresponding to the spectrum above c (the latter 
is not bounded above). By Theorem l2.31 we know that using such a basis will not create any 
pollution in the gap of A. 

We emphasize that the Wannier basis does not depend on the decaying potential W, and 
can be precalculated once and for all for a given ^ and a given Vper- Another huge advantage 
is that since Wk decays fast, it will be localized over a certain number of unit cells of .if. As 
W represents a localized defect in the lattice, keeping only the Wannier functions Wk{- — a) 
with a S .5f ni?(0, R) for some radius R > should already yield a very good approximation 
to the spectrum in the gap (we assume that the defect is localized in a neighborhood of 0). 
This approximation can be improved by enlarging progressively the radius R. 

Of course in practice exponentially localized Wannier functions are not simple to calcu- 
late. But some authors have defined the concept of maximally localized Wannier functions 
[24] and proposed efficient methods to find these functions numerically. 

The efficiency of the computation of the eigenvalues of A in the gap using a Wannier 
basis (compared to that of the so-called super-cell method) were illustrated for a nonlinear 
model in [6j. 



2.3.2. Dirac operators in upper/lower spinor basis 

The Dirac operator is a differential operator of order 1 acting on L^(R^,C'*), defined as [38l 

US] 

3 

D° = -ic^akd^^ +mc^l3 := ca-p + mc"^!]. (2.19) 

k=l 

Here ai,a2, as and f3 are the so-called Pauli 4x4 matrices [38^ which are chosen to ensure 
that 

{D° f ^-c^A + m^c\ 
The usual representation in 2 x 2 blocks is given by 

^-i'oX)' "-C/o) '*-'-^.»)' 

where the Pauli matrices are defined as 

-(■;)■ 

In the whole paper we use the common notation p = — «V. 

The operator is self-adjoint on H^{M.^, C^) and its spectrum is symmetric with respect 
to zero: cr(D°) = ( — cxD, —mc^] U [toc^, oo). An important problem is to compute eigenvalues 
of operators of the form 

^D^^ + V 

in the gap (— mc^, mc^), where y is a multiplication operator by a real function x i— > y{x). 
Loosely speaking, positive eigenvalues correspond to bound states of a relativistic quantum 
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electron in the external field whereas negative eigenvalues correspond to bound states of 
a positron, the anti-particle of the electron. In practice, spectral pollution is an important 
problem [131 [HI (HI [34] which is dealt with in Quantum Physics and Chemistry by means 
of several different methods, the most widely used being the so-called kinetic balance which 
we will study later in Section [3.2.11 We refer to [3] for a recent numerical study based on 
the so-called second-order method for the radial Dirac operator. 

We now present a heuristic argument which can be made mathematically rigorous in 
many cases [38l [15]. First we write the equation satisfied by an eigenvector {(p, x) of D'^ + V 
with eigenvalue mc^ -t- A e (— mc^, mc^) as follows: 



(mc2 + V)(p + CO- ■ (-iV)x = (mc2 + X)ip, 
{—mc^ V)x + C(T • {—iV)ip — [mc^ + X)x 



where we recall that a = (ci, 0-2, fa) are the Pauli matrices defined in (|2.20p . Hence one 
deduces that (when it makes sense) 

(2.22) 

If V and A stay bounded, we infer that, at least formally. 



X 



TT^cr • {~iV)ip 
2mc ^ 



(2.23) 



Hence we see that in the nonrelativistic limit c — s- 00, the eigenvectors of A associated with 

[•A ■ ■ 

a positive eigenvalue converge to a vector of the form _ . Reintroducing the asymptotic 



formula (|2.23p of x in the first equation of (|2.2ip , one gets that if is an eigenvector of the 
nonrelativistic operator —A/ {2m) + V in L^(M^,C^). 

For this reason, it is very natural to consider a splitting of the Hilbert space L^(R^,C'*) 
into upper and lower spinor and we introduce the following orthogonal projector 



(2.24) 



This splitting is the choice of most of the methods we are aware of in Quantum Physics and 
Chemistry. Applying Theorem 12. 1[ we can characterize the spurious spectrum associated 
with this splitting. For simplicity we take m = c = 1 in the following. 

Theorem 2.4 (Pollution in upper/lower spinor basis for Dirac operators). Assume 
that the real function V satisfies the following assumptions: 

(i) there exist {Rk}kLi ^ o,''^d a positive number r < inffc^^ \Rk — Re\/2 such that 

/3 

max sup \x- Rk\\V{x)\ < ^] (2.25) 

(u) one ha^ 

Vtm\u'<B(R,.r) e L'^iR') + LT(M') for some 3 < p < oo. (2.26) 



=We use the notation of ^O,: X + = {f e X + L°° \ \/€ > 0, 3fc e X such that ||/ - /el^oo < e}. 
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Let V he as in ()2.24|) . Then one has 

Sp^zFTV^ = {Conv (Ess(l + V)) U Conv (Ess( - 1 + V)) } n [-1, 1] (2.27) 

where Ess(W^) denotes the essential range of the Junction W , i.e. 

Ess(VF) = {A e R I |T^~\[A - e, a + e])| 7^ Ve > O}. 

Remark 2.6. It is known that the operator + V is essentially self-adjoint on C^(K.^, C^) 
when (j2.25p and (I2.26P hold, and that its domain is simply the domain _ff^(]R'^,C^) of the 
free Dirac operator. When \/3/2 is replaced by 1 in (|2.25p . the operator + V still has 
a distinguished self-adjoint extension 138] whose associated domain satisfies iJ^(R'^,C^) C 
D(D^ -\- V) G -ff^/^(R'^, C*). Furthermore this domain is not stable by the projector V on 
the upper spinor (a characterization of this domain was given in llGp. The generalization 
to this case is possible but it is outside the scope of this paper. 

Remark 2.7. By Theorem [g^j we see that Spu{D°,V) = but Spu{D° + V,V) for 
all smooth potentials V even if V is -compact. Hence spectral pollution is in general 
not stable under relatively compact perturbations (but it is obviously stable under compact 
perturbations as we have already mentioned in Remark l 1 . S^) . 

Our assumptions on V cover the case of the Coulomb potential, V{x) — k\x\^^ when \ k\ < 
^/3/2. In our units, this corresponds to nuclei which have less than 118 protons, which covers 
all existing atoms. On the other hand, a typical example for which V G Lp(R^) n L°°(R'^) 
is the case of smeared nuclei V — p*l/\x\ where p is a (sufficiently smooth) distribution of 
charge for the nuclei. We now give the proof of Theorem 12.41 

Proof. Under assumptions (z) and (m), it is known that D'^ + V is self-adjoint on i/^(R'^, C^) 
(which is stable under the action of V) and that cross(-D'^ -\-V) — (— oo, —1] U [1, oo). We will 
simply apply Theorem 1 2 . 1 1 wit h C = H^{M.'^,C^). We have in the decomposition of L^(R'^) 
associated with V, 



Hence ViD° + V)V = 1 + F and (1 - V){D^ + V){1 - V) = -l + V, both seen as operators 
acting on L^{R^,C^). It is clear that D{D" + V) n VL^iR^,C'^) ~ H^R^,C^) is dense in 
the domain of the multiplication operator by V{x) 




D{V) = {f e L^{R^,C^) I Vf e L=^(R^C2)}, 



for the associated norm 




Also the spectrum of V is the essential range of V. Note under our assumptions on V we 
have that G Ess(y). The rest follows from Theorem 12. II □ 
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2.3.3. Dirac operators in dual basis 

In this section we study a generahzation of the decomposition into upper and lower spinors, 
which was introduced by Shabaev et al [32]. For any fixed e, we consider the unitary operator 

which is just a dilation of the sign of Z?" (note that (U^)* = U^). Next we define the following 
orthogonal projector 

Ve ■■= UeVU, (2.29) 

where V is the projector on the upper spinors as defined in (|2.24[) . As for e = we have 
Uq = 1, we deduce that Vq = V . However, as we will see below, the limit e ^ seems to 
be rather singular from the point of view of spectral pollution. We note that any vector in 
PeL^(]R'^, C*) may be written in the following simple form 



with LP e H^i^?,^'^ 



ea ■ {—iV)Lp ^ 

Hence for e ^ 1, the above choice just appears as a kind of correction to the simple decom- 
position into upper and lower spinors. Also we notice that VeH^{M.'^,C^) C ff^(K^,C'*) for 
every e since Ue is a multiplication operator in Fourier space and V stabilizes H^{MP,C'^). 

In [32], the projector Ve is considered with e = l/(2mc) as suggested by Equation ()2.23p . 
However here we will for convenience let e free. The method was called "rfMaf in |32j since 
contrarily to the ones that we will study later on (the kinetic and atomic balance methods), 
the two subspaces VeL'^lR^ ,C^) and (1 - 7'<;)i^(R^, C) play a symmetric role. For this 
reason, the dual method was suspected to avoid pollution in the whole gap and not only in 
the upper part. Our main result is the following (let us recall that m = c = 1): 

Theorem 2.5 (Pollution in dual basis). Assume that the real function V satisfies the 
following assumptions: 

(i) there exist {Rk}kLi C K"^ and a positive number r < inf^^^ \Rk — Ri\/2 such that 

/3 

max sup \x~ Rk\\V{x)\ (2.30) 

''='^--^ \x~B.k\<r 2 

(ii) one has 

^lR3\ufB(i?fc,r) e £P(R^)nL°°(M^) for someS <p<oo. (2.31) 
Let < e < 1 and Ve as defined in (|2.29p . Then one ha^ 



Spu(Z3" + y,7',) = 



2 

— 1 , min ^ — - + 1 + sup V , 1 



4 , ^-l + infyj. , 1 



(2.32) 

Our result shows that contrarily to the decomposition into upper and lower spinors 
studied in the previous section, the use of Ve indeed allows to avoid spectral pollution under 
the condition that is a bounded potential and that e is small enough: 

2 



"^Recall that [a, fe] = if 6 < a. 
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This mathematically justifies a claim of |32j . However we see that for Coulomb potentials, 
we will again get pollution in the whole gap, independently of the choice of e. Also for large 
but bounded potentials (like the ones approximating a Coulomb potential) , one might need 
to take e so small that this could give rise to a numerical instability. 



Proof. We will again apply Theorem O We choose C = U^Cg° {M.^ , C^) . Note that C is a 
core for D° + V (its domain is simply H^(R^,C^)) and that VeC C ns>oi?^(M3, C") since 
Ve and Ue commute with the operator p = — iV. An easy computation yields 



2 



-y)^^^---A, (2.33) 



and 



C/,(D" + F)|(i_p,)cC/e =^ -1 + , „ V 



'"•^ ^-^2 + ^) /"-P :=A,. (2.34) 



Strictly speaking these operators should be defined on VUeC and (1 — V)UeC but we have 
made the identification VU^C ~ (1 - 'P)UeC ~ C^(K3,C2). Let us remark that for e > 
the term K := {1 + e^|pp)~^/^V^(l + e^|pp)~^/^ is indeed compact under our assumptions 
on V, hence it docs not contribute to the polluted spectrum. On the other hand for e = it 
is the only term yielding pollution as we have seen before. 

Theorem 12.51 is then a consequence of Theorem 12.11 and of the following 



Lemma 2.3 (Properties of {D'^ + V)^-p^c (Z)" + V')|(i_-p^)c). The operators Ai and 
A2 defined in (|2.33p and (j2.34p are self-adjoint on the domain 

P := e L\R\C^) I V{a-p){l + e^\p\^)-^^^ip G L''{W^,C^)} . 

They are both essentially self-adjoint on C^(R"^,C^). Moreover, we have 

Conv Ess ( - - 1 + y ) C Conv cToss (^1) ^ 



e 



C 



,2 1 f 2 

min <j 1 , 1 + inf V > , max < 1 , 1 + sup V 



an I 



d 



Conv Ess ( - - + 1 + ) C Conv CTcss (^2) C 



e 



C 



2 1 r 2 

-1 , h 1 + inf \/ ^ , max -^-1, hl + sup^ 



Proof. The operator K — {1 + e^|pp)~^/^y(l + e^|pp)^^/^ being compact, it suffices to 
prove the statement for £e(— 2/e + 2 + V)£i, where we have introduced the notation := 
ea ■ p{l + e^lpp)^^/^. The argument is exactly similar for L^{2/e — 2 + V)Lc. We denote 
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W :— —2/e + 2 + V and we introduce A = C^WC^ which is a symmetric operator defined 
on T). We also note that T> is dense in L^. 

Let / e D{A*), i.e. such that \{f, CeWCe(p)\ < C\\ipl ^tp G V. We introduce x ■= 
lyM ^B(Rk r)j ^ localizing function around the singularities of and we recall that V is 
bounded away from the i?fc's. Hence we also have KfjC^x^^ef)] ^ C" \\ip\\ for all ip £ V. 
Then we notice that under our assumptions on V, we have Wx € L^, hence g := x^^ef S 
L-^ and g G L°°. In Fourier space the property \JgC^(p\ < C \\(p\\ for all (p in a dense 
subspace of means that ecr • p{l + e^|pp)"^/^5(p) G i^, hence 5 G L^(M^ \ B(0, 1)). As 
by construction g G we finally deduce that g G L^, hence WC^f G L^. We have proven 
that D{A*) C I?, hence A is self-adjoint on V. The essential self-adjointness is easily verified. 

The next step is to identify the essential spectrum of A. We consider a smooth normalized 
function C G C|^(R^,]R) and we introduce ipi = {1 + a ■ p/|p|)(C,0). We notice that (pi G 
_ff'*(R^, C^) for all s > 0. Then we let (pnix) := n'^^^ipi{n{x—xo)) andnotetha.t {<7-p/\p\)(pn = 
(pn- We take for xq G 



1^ some fixed Lebesgue point of V, i.e. such that 
1 



lim 



r^O \B{xo,r)\ JBixo,r) 



\V{x) - V{xo)\dx = 0. 



(2.35) 



First we notice that 



lim 



ecr • p 



-l] ipn 



= lim 



1 <y5r: 



= 



as is seen by Fourier transform and Lebesgue's dominated convergence theorem. Therefore, 



lim 

n — >oo 



w 



ecr • p 



1 Cn 



L2 



since we have W E L'^ + L°°. On the other hand we have lim„ 



.WiW-WixoJKn 



(2.36) 



= 0. 



. \\{A - W{xo)Kn\\L2 = 0. This proves 

^ CSS 

(A) as seen by taking (Pn{x) = 



Using this to estimate cross terms we obtain lim„^( 
that Ess(W^) C cress (^)- Let us remark that G 
n^'^^^ipi{x/n). 

The last step is to show that CTcss(^) Q [min{0, inf(P^)}, max{0, sup(P^)}]. When 
sup(W) < 00, we estimate A < sup{W)C^. If < 0, then we just get A<0, hence a{A) C 
(—00,0]. If < sup(VK) < 00, we can estimate £^ < 1 and get <t{A) C (— 00, sup(lF)]. 
Repeating the argument for the lower bound, this ends the proof of Lemma [2.31 □ 



2.3.4. Dirac operators in free basis 

In this section, we prove that a way to avoid pollution in the whole gap is to take a basis 
associated with the spectral decomposition of the free Dirac operator, i.e. choosing as pro- 
jector :— X{o.oo){D^)- As we will see this choice does not rely on the size of V like in 
the previous section. Its main disadvantage compared to the dual method making use of Ve, 
is that constructing a basis preserving the decomposition induced by requires a Fourier 
transform, which might increase the computational cost dramatically. First we treat the case 
of a 'smooth' enough potential. 

Theorem 2.6 (No pollution in free basis - nonsingular case). Assume that V is a 
real function such that 

V G lp(r3) + (l'^(r3) n w^'''(m.^)^ + Lf(m.^) 
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for some 6 < p < oo, some 3 < r < 6 and some 2 < q < oo. Then one has 

Spn{D" + V,Pl) 
Remark 2.8. We have used the notation 

l''(m3) n w^'iiR^) = {V e ^'■(M^) I e L'^{w^)}. 

Remark 2.9. A physical situation for which the potential V satisfies the assumptions of 
the theorem is V ^ p * with p G ^^(R^) n L'^{W'). 

Proof. Under the above assumptions on the potential V , it is easily seen that the operator 
D° + y is self-adjoint with domain iJi(R3,C''), the same as , and that a^ssiD^ + V) = 
c{D^) — (— oo, — 1] U [1, oo) (these claims are indeed a consequence of the calculation below). 
Hence (—1, 1) only contains eigenvalues of finite multiplicity of + V and we may find a 
c G (—1, f) \ a{D^ + V). In the following we shall assume for simplicity that c = 0. The 
argument is very similar if G a{D^ + V). We will denote 11 — X[q,od){D^ + ^) a-nd prove 
that {Pi - IV)\D^ + is compact. This will end the proof, by Theorem [2:21 

As a{D° + V), we have that + V\>e for some e > 0. Also, we have 

for e > small enough. Taking the square root of the above inequality, this proves that 
|£)0|-i/2|£)0 _|_ ^|i/2 g^jj^j j^g inverse are both bounded operators. 

Next we use the resolvent formula together with Cauchy's formula like in [19] to infer 

(P" - n)\D° + V\^/^ =-— r ( — ^ \D" + V\'/^drj 

1 1 ^\D° + V\^/^ , 



27r D" +iri + V + iri 

Let us now write ¥ = + + K," with V{' G LP{R^) for 6 < p < 00, G L 
and VVi" G for 3 < r < 6, 2 < 9 < 00, and jVi'W^^^j^s) as n 00. We write 



(p° - n)|D° + ^ K{v{') + K{v:^) + k{v;') 



with 



1 r 

K{W) := - / 



1 rr,\D° + V\^, 



D° + iT] + V + 

and estimate each term in an appropriate trace norm. We denote by 6p the usual Schatten 
class [33l[30] of operators A having a finite p-trace, \\A\\q = Tr{\A\Py^P < 00. Let us recall 
the Kato-Seiler-Simon inequality (see [31] and Thm 4.1 in [33] ) 

Vp > 2, ||/(-zV)5(x)||e^ < (27r)-3/f ||/||^,(r3) ||<7|Il.(e3) ■ (2.37) 

The term K(V{^) is treated as follows: 

1 r *7_ 

2^ J -00 (£2+772)1/4 



where we have used that 

lD" + X/|i/2 ^ 









- y + i?7 



(£2+^2)1/4- 
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By ([237]) we have 



dk 



|fc|2 + 772)p/2 



1/p 



< 



c 



1 + 77 



||T/n| 



(2.38) 

, hence this term is a 



Since 6 < p < 00, this finally proves that \\K{V{')\\q^ < C 
compact operator for any n. 

The term involving is more complicated to handle. First we use the formula [121 HO] 



1 



-V, 



1 



D° + ir] ^ D° +iri 



di] = 



1 



DO + IT] 
1 

1 



1 



2 J 



+ iri 



dr/ 



+ iri 



dr] 



(DO + ir;)- 



r(«-Vl/2" 



+ iri 



— dr]. 



Iterating the resolvent formula we arrive at 



(DO + irj) 



^ (a . WV-)^^^dv ] |D°rV2|^o ^ ^|i/2 



D" + Z77 



1 



1 1 |D0 + 1/|1/2 

^2-7^7^^ ' n , . dv. (2.39) 



27r + iri ^ D" + ii] + V + ir] 

The first term can be estimated as before by (recall that |D0|^i/2|£)0 y|i/2 jg bounded) 



2^ 



(DO + ir;) 



1 |D0|1/2 



DO + i77 



< CIIVK 



2 IIL9(R3) 



dry 



-00 1 + r; 



1+3^ 



which is convergent since g > 2 by assumption. 

The next step is to expand the last term of (|2.39p using again the resolvent expansion: 



1 1 |D0 + 1/|1/2 



+ ir] ^ + ir] D° + V + ir] 

fc-i 



ir] 



DO 



(-1)' 



1 



+ ir] J DO + 11] I 

|D° + F|i/2dr/. (2.40) 



_^D" + ir] ^ \DO + ir] J D^ + V + ir] 



By (I2.38P we see that the last term belongs to 6kr when k is chosen large enough such that 
fc(l — 3/r) > 1/2 (which is possible since r > 3). 

We now have to prove that the other terms corresponding to j = 2...fc — 1 in (|2.40p are 
also compact. We will only consider the term j = 2, the others being handled similarly. 
Writing V = + V2 + V3" the terms containing and are treated using previous 
ideas. For the term which only contains V2 , the idea is, as done previously in [19], to insert 



po + po = 1 as follows 



00 pO I pO pO I pO pO I pO 

^+ + - T/« _±_L_T/" + - I nO 1 1/2^„ 
DO + zry ^2 ^0 + >'2 ^0 + 1^ I ^'^^ 
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The next step is to expand and note that, by the residuum formula, the ones which contains 
only or only vanish. Hence we only have to treat terms which contain two different 
P^. We will consider for instance 

P° . Pi Pi 



\rn '+ yn + |n0|l/2n'„ 



/oo pU pU pU 

73ir^[^°,^2"]7«rf-^2"75jrf-l^°l'/'rf^- (2.4i) 
+ 17] + irj D" + 

Now we have using again a Cauchy formula for Pi 

■ poo |r)0|-l/2 1 

The Kato-Seiler-Simon inequality (|2.37p yields as before 



< CIIVKT 



Inserting this in ()2.4H) and using that VJ^ G L*", we see that the corresponding operator is 
compact. 

Eventually we have by a trivial estimate 

{Pi - n)|i?" + _ K{vn _ K{v^) ^ \\K{v^)\\ < c Wv^l^^^,^ ^„^oo o. 

As (P° — n)|Z3° + V\^^^ is a limit in the operator norm of a sequence of compact operators, 
it must be compact. □ 

We treat separately the case of a Coulomb- type singularity, for which (P° — n)|£''^ + y|^/^ 
is not compact, hence we cannot use Theorem 12.21 directlv. 

Theorem 2.7 (No pollution in free basis - Coulomb case). Let \k\ < \/3/2. Then 

Spu (^7^o + ^,P^^n(-i,i) = 0. 

Proof. The operators Pl{D" + k\x\^^)PI and P"{D° + k\x\^^)P1 arc known to have a 
self-adjoint Friedrichs extension as soon as \k\ < 2/{tt/2 + 2/tt), see |17j. Furthermore one 
has (Tess(-D° -I- k|x|~^)|po^2 = [l,oo) and tToss(£'° + K|a;|~^)|po ^2 — (— 00, —1], see Theorem 

2 in [13. As ^/3/2 < 2/{tt/2 + 2/tt), the result immediately follows from Theorem [2T] and 
Remark [1?T] □ 



3. Balanced basis 

In Section [2] we have studied and characterized spectral pollution in the case of a spitting 
Sj = Pjo® (1 — P)i] of the main Hilbert space. In particular for the case of the Dirac operator 
+ y we have seen that the simple decomposition into upper and lower spinors may yield 
to pollution as soon as 7^ 0. In this section we study an abstract theory (inspired of 
methods used in Physics and Chemistry) in which one tries to avoid pollution by imposing 
a relation between the vectors of the basis in P9j and in {1 — P)f), modelled by one operator 
L : PS) (1 — P)Sj. We call such basis a balanced basis. 
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3.1. General results 

Consider an orthogonal projection P : 9j ^ Sj. Let L : D{L) C PSj ^ (1 — P)Sj be a 
(possibly unbounded) operator which we call balanced operator. We assume that 

• L is an injection: if Lx — for x G D{L)^ then a; = 0; 

• D{L) ® LD{L) is a core for A. 

Definition 3.1 (Spurious eigenvalues in balanced basis). We say that A G M is o 
(P, L)-spurious eigenvalue of the operator A if there exist a sequence of finite dimensional 
spaces {V^}n>i C D{L) with C V^j^i for any n, such that 

T 1-D{A) 

(1) \Jn>i{V,t ® LV^) = D{A); 

(2) £m dist {\<y{\v+®LV^))) = 0; 

(3) \ia{A). 

We denote by Spu(yl, P, L) the set of {P, L)- spurious eigenvalues of the operator A. 

Remark 3.1. Another possible definition would be to only ask that for all n, V^^ contains 
LVj^ . This would actually also correspond to some methods used by chemists (like the so- 
called unrestricted kinetic balance mV - The study of these methods is similar but simpler 
than the one given by Definition \3.1\ 

Contrarily to the previous section, we will not characterize completely (P, i)-spurious 
eigenvalues. We will only give some necessary or sufficient conditions which will be enough 
for the examples we are interested in and which we study in the next section. We will assume 
as in the previous section that PAP (resp. (1 — P)A{\ — P)) is essentially self-adjoint on 
D{L) (resp. on LD{L)) with closure denoted as A|p^ (resp. A|(i_p)f|). 



3.1.1. Sufficient conditions 

We start by exhibiting a very simple part of the polluted spectrum. For any fixed ^ 
X G D{L), we consider the 2x2 matrix M{x) of A restricted to the 2-dimensional space 
xC ® LxC, and we denote by fJ,i{x) < H2{x) its eigenvalues. Note that jii is homogeneous 
for 4 = 1,2, Hi{\x) = ^ii{x). 

Theorem 3.1 (Pollution in balanced basis - sufficient condition). Let A, P, L as 

before and define nii, Mi G M U {±cxd}, i = 1, 2, as follows: 

mi:— inf liminf /ii(a;^). Mi:— sup limsup/xi(x^), (3.1) 

{xi}CDiL)\{0}, {xt}CD{L)\{0}. »^°° 

Lxt\\Lxif^O Lxt\\Lxt\\''^0 



m2 :— inf liminf /X2(a;^), M2 := sup limsup/i2(a;^). (3.2) 

{xt}CDiL)\{0}, {:i;+}cr>(L)\{o}, 

L<\\Lxt\\~'^0 Lxt\\Lx+\\-'^0 



Then we have: 



[mi, Ml] U [m2, M2] C Spu(^, P, L) U acssiA). (3.3) 



26 



Mathieu LEWIN & Eric SERE 



We supplement the above result by the following 

Remark 3.2. The two diagonal elements of the matrix A(x^) being {Ax^,x^) \\x^\\ ^ and 
{ALx^,Lx^) \\Lx^\\ ^, it is clear that we have 

m2>mi> max(inf(Tcss(^|(i-P)jj), inf <Toss(A|p.(j)) , 

Ml < M2 < min (sup (Toss (A|(i_p).(j) , SUp(Tcss(^|Pf,)) , 

which is compatible with Theorem \2.1[ since we must of course have Spu(A, P, L) C 
Spu(A,P). 

Proof. We will use the following 

Lemma 3.1. Assume that A, P and L are as above. Let {Vn} C D{L) be a sequence of K- 
dimensional spaces with orthonormal basis (x^, x^). Let be an orthonormal 

basis of LVn C (1 — P)Sj. We assume that xj^ and y'^ weakly for every k — 1..K , as 
n 00. If X E R is such that lim„^oo dist (A , aiA^y )) = 0, then A G Spu(A, P, L) U 
a{A). 

The proof of Lemma 13.11 will be omitted, it is very similar to that of Lemma 11.21 We 
notice that the two sets 

:= |/i e MU {±00} : 3{x+}(lD{L), a;+ ||a;+||"^ ^ 0, 

L4||L4f' -0, (3.4) 

are closed convex sets, for i — 1,2. Indeed, assume for instance that Ai, A2 G Ki and let be 
{xn\ and {j/n} such that /ii(a;„) ^ Ai and tiiiyn) — > A2. By the homogeneity of /ii we may 
assume that ||x„|| = ||y„|| = 1 for all n. Also, extracting a subsequence from we may 

always assume that 

lim {xn,yn) ^ lim ( 77-^' ) = 0- 



n — >-oo 



^ \Lxn\ ' \Ly„ 

Fix some A G (Ai, A2) and consider as usual Zn(9) = cos 6 x„ + svaO yn- By continuity 
of the first eigenvalue of the 2x2 matrix of A in the space spanned by ZniO) and Lzn{9), 
we know that there exists (for n large enough) a 0„ G (0, 27r) such that /Zi(0„) = A. Note 
that ||z„(6l„)|| = 1 + 0(1). Writing Lz„(6'„) \\Lzn{On)\\~^ = a„ia;„ ||La;„||~^ + (3nLyn ||i?;„||~"^ 
we see that both a„ and /3„ are bounded and satisfy 0^ + 0^ 1, hence ||Lz„(6'„)|| — > 1. 
It is then clear that Zn{0n)\\zn{6n)\\~^ ^0 and that Lz„(0„)||Lz„(0„)|| ^ 0. Therefore 
A = lim„^oo fJ'2{zn{&n)) G Ki- The argument is the same for K2. As Lemma l3. II tells us that 
Ki U K2 C Spu{A, P, L) U a{A), this ends the proof of Theorem O □ 



3.1.2. Necessary conditions 

Let us emphasize that, contrarily to P-spurious eigenvalues, for (P, L)-spurious eigenvalues 
the two spaces Pio and (1 — P)f) do not play anymore a symmetric role due to the in- 
troduction of the operator L. For this reason we shall concentrate on pollution occurring 
in the upper part of the spectrum and we will not give necessary conditions for the lower 
parlE. Loosely speaking, obtaining an information on the lower part would need to study 

•^As wc have mentionned before we always assume for simplicity that inf (Tess(^| (i-p)fj ) < inf o"css(A|pjj ), 
i.e. that 1 — P is responsible from the pollution occuring in the lower part of the spectrum. 
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the operator L ^ . In the applications of the next section, we will simply compute the lower 
polluted spectrum explicitely using Theorem 13. II Let us introduce 



d := supcr(A(i_p)jj). (3.5) 

and assume that d < oo. In the sequel we will only study (P, L)-spurious eigenval- 
ues in (c, oo). Note that due to Theorem 12. li it would be more natural to let instead 
d := sup (Toss (^(i-pjij) but this will actually not change anything for the examples we 
want to treat: in the Dirac case I?" + V and for P — V, the orthogonal projector on the 
upper spinor defined in (|2.24[) . the spectrum of (D° + T^)|(i_p)l2(]j3 ^4) = —1 + V is only 
composed of essential spectrum. We do not know how to handle the case of an operator 
A|(i_p)i5 which has a nonempty discrete spectrum above its essential spectrum. Our main 
result is the following 

Theorem 3.2 (Pollution in balanced basis - necessary conditions). Let A, P, L as 

before. We recall that the real number d < oo was defined in p.Sp . 

(i) Let us define 

m'n = inf f'2ix^) (3.6) 

x+6D(L)\{0} 

and assume that rrij > d. Then we have 

SpvL{A,P,L) n (rf,TO2) = 0- 

(ii) Let us define 

777.2= iiif lim inf /i2 (x^ ) (3.7) 

{^i}cr'(L)\{o}, 

x+^O, ||a;+|| = l 

and assume that m'2 > d. We also assume that the following additional continuity property 
holds for some real number b > d: 



{x+} C DiL) 
limsup/i2(2^^) < b 



{Axt,x+)^0. (3.8) 



Then we have 



Spu(A, P, L) n {d, min(7772, 6)) = 



Remark 3.3. The property ()3.8|) is a kind of compactness property at of the set {x^ G 
D[L) I ^2(2;^) < 6} for the quadratic-form norm of the operator A\ps). 

Remark 3.4. Note that (|3.8p holds true for b — +00 > d when A^p^^ is a bounded operator. 



Theorem l3.2l has many similarities with the characterization of eigenvalues in a gap which 
was proved by Dolbeault, Esteban and Sere in [12j (where our number d — sup (T(^(i_p)ij) 
was denoted by 'a'). In particular the reader should compare the assumptions d < m'2 and 
d < m'2 with (iii) at the bottom of p. 209 in [12]. The proof indeed uses many ideas of 
[12] . Note that [12] was itself inspired by an important Physics paper of Talman [37] who 
introduced a minimax principle for the Dirac equation in order to avoid spectral pollution. 
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Proof. Assume that A G Spu{A, P, L) n {d,oo). We consider a Weyl sequence {a;„} like in 
Lemma ri-H i.e. such that 

P^+^^^+(A-A„)x„ =0 (3.9) 

for some x„ ^ with ||x„|| = 1 and some A„ — s- A. We write Xn ~ where x^ = Ly^ 

for some j/+ G . Now, hke in [12j we consider the following functional defined on LV^: 

Q{x~) := {A{x^ + + x^) - A„ + x^\\^ . 

Using the equation Py+^j^y+{A — A„)a;„ — 0, we deduce that 

Vx^ £ LV+, Q{x^) = {{A - \n){x~ - x~),x^ - x:^). 

By definition of d we obtain 

\/x-eLV,t, Q{x-)<{d-\n)\\x- -x-f . (3.10) 

Consider the 2x2 matrix M{xf^) of A restricted to x^ © ix+ and recall that fi2{Xn) is 
by definition its second eigenvalue, hence 

+ {Aix+ +eLx+),x+ +eLx+) 

M2«) = sup- 



(the sup is not necessarily attained). There exists 0„ £ R such that for n large enough 

j— — 12 >M2«)-l/n- (3.11) 

Inserting = O^Lx^ in (j3.10p we obtain for n large enough, 

(M2(a:+) - A„ - 1/n) f + 0^ f ) + (A„ - d) - a." f < 0. (3.12) 

Let us assume we are in case (i) for which mi^ > d. Using the obvious estimate /i2(a;^) > 

we see that if A € (c?, mj), then for n large enough we must have x^ = 9nLx^ = 6nLx^ — 
x~ =0,thusa;„ = which is a contradiction with ||a;„|| = 1. Hence Spu(A, P, L)n(d, rrij) = 0. 

Let us now treat case (ii) for which we assume TO2 > d and that p.Sp holds for some b > d. 
Let A e Spu(A, P, L)n(d, min(6, mj)). From (|3.12p we see that necessarily H2{Xn) < A„ + l/n 
(except if a;„ = which is a contradiction) . Therefore we have lim sup^^^^ /12 (a;^ ) < b. 

Assume first that a;^ ^ strongly. Using our assumption (|3.8|) . we deduce that 
lim„^oo {Ax^,x^) = 0. Next we argue like in the 3rd step of the proof of Theorem 12.11 
First, taking the scalar product of (|3.9|) with we deduce that lim„^oo (Ax^^x^) — 0. 
Taking then the scalar product with x^ we deduce that 

lim {{A - A„)a;,7,a:;,7) = 0. 

n— i-oo 

As {{A — Xn)x~ ,x~) < (d — A + 0(1)) \\x~ II and d — A < we deduce that x:^ ^ which is 
a contradiction with ||a;„|| = 1. 

Hence we must have x^ ^ 0, which implies that x^ W^nW ^ 0, up to a subsequence. 
Therefore we have liminf„^oo /i2(a;^) = lim inf „^oo ^^2 (a;^ ||a;^|| ^) > m'2 by definition of 
7712 . Inserting this information in (j3.12p . we again arrive at a contradiction, similarly as 
before. This ends the proof of Theorem 13.21 □ 
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3.2. Application to Dirac operator 

In this section, we consider the Dirac operator A — + V for a potential satisfying the 
assumptions (|2.25p and (|2.26|) of Theorem [2^ and 

sup(y) < 2 (3.13) 

We wiU indeed for simphcity concentrate ourselves on the case for which either V is bounded, 
or F is a purely attractive Coulomb potential, V{x) — — I, < K < V3/2. The general- 
ization to potentials having several singularities is rather straightforward. 

Like in Section [2.3.21 we start by choosing P = V, the projector on the upper spinors 
as defined in (|2.24p . As already noticed in Section [2.3.21 we then have VAV — 1 + V and 
(1 — V)A{1 — V) = — 1 + y on the appropriate domain. This shows that the number d 
introduced in the previous sections is d = — 1 + supF < 1 by (|3.13p . 

Wc will now study different balanced operators L which we have found in the Quantum 
Chemistry litterature. Note that we can always see L as an operator defined on 2-spinors 
D{L) C L'^{R^,C'^) with values in the same Hilbert space i^(R'^,C^), which we will do in 
the rest of the paper. 

We will describe the polluted spectrum Spu(-D*' + VjV, L) using the results presented in 
the previous sections. We note that the number iJL2{'f) is the largest solution to the following 
equation [12] 

«^ + ^> + {{^. + l-v)L^,L^) ^ ^ 

where the denominator of the second term does not vanish when ijl2['^) > d = sup(y) — 1. 
Note the term on the left is decreasing with respect to n, whereas the term on the right is 
increasing with respect to /z. Hence we have IJ-2{'P) > 1 if and only if 

where the denominator of the second term does not vanish due to p.l3p . Note that (I3.15P 
takes the form of a Hardy- type inequality similar to those which were found in [l2l [11] . In 
the following wc will have to study this kind of inequalities for sequences which converge 
weakly to 0. The Hardy inequalities of [12|[Tl] will indeed be an important tool as we will 
see below. 

Concerning the choice of the operator L, several possibilities exist, although the main 
method is without any doubt the so-called kinetic balance which we will study in the next 
section. All the methods from Quantum Chemistry or Physics are based on the following 
formula for an eigenfunction {(p, x) with eigenvalue mc^ -t- A (we reintroduce the speed of 
light c and the mass m for convenience) and which we have already formally derived before 
in Section [2X1 

This equation suggests that for an eigenvector to be represented correctly, the basis of the 
lower spinor should contain c{2mc^ -I- A — V)~^a ■ (— iV) applied to the elements of the basis 
for the upper spinor. However we cannot choose in principle L = c{2mc^ -t- A — V)~^a- (— iV) 
because A is simply unknown. For this reason, one often takes the first order approximation 
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in the nonrelativistic limit which is nothing but 

Lkb = 77— o- • (-«V). 
2mc 

The choice of this balanced operator is (by far) the most widespread method in Quantum 
Physics and Chemistry. It will be studied in details in Section fS. 2. II 

It seems a well-known fact in Quantum Chemistry and Physics [111 [27] that the kinetic 
balance method consisting in choosing L = Lxb is not well-behaved for pointwise nuclei. 
The reason is that the behaviour at zero of c{2mc^ + A — V)~^a ■ (— «V) is not properly 
captured by a ■ (— iV), if V{x) = —k\x\~^. Indeed we will prove that the kinetic balance 
method allows to avoid pollution in the upper part of the spectrum for 'regular' potentials, 
but not for Coulomb potentials, which justifies the aforementioned intuition. 

To better capture the behaviour at zero, we study another method in Section [3 . 2 . 21 which 
we call atomic balanc^ and which consists in choosing 

= 2^^^^" ■ 

Although this operator does not depend on A, it will be shown to completely avoid pollution 
in the upper part of the spectrum, even for Coulomb potentials. It is very likely that any 
other reasonable choice with the same behaviour at zero would have the same effect but we 
have not studied this question more deeply. 

In the following we again work in units for which m — c — 1. 

3.2.1. Kinetic Balance 

The most common method is the so-called kinetic balance [T3|, [181 [221 [34] . It consists in 
choosing as balanced operator 



j^B ^-ia-V (3.17) 



We can for instance define Lkb on the domain D{Lkb) — C^{M? , C^), in which case Lkb 
satisfies all the assumptions of Section [3] Our main result is the following 

Theorem 3.3 (Kinetic Balance), (i) Bounded potential. Assume that V E LP{M.^) 
for some p > 3, that lim|j.|^oo ^{^) = 0; '"^'^ ^^'J^ 

- 1 + sup(V^) < 1 + inf(l/). (3.18) 

Then we have 



Spn{D° + V,V,Lkb) = [-1,-1 + supF]. 
(ii) Coulomb potential. Assume that < k < \/3/2. Then we have 



Spul D" - -l^,V, Lkb] = [-1,1]. (3.19) 



Remark 3.5. The conclusion ([?1^ also holds if V is such that V G LP{R^) n L°°(R^ \ 
B{xq, r)) for some p > 3 and 



< V{x) < on B{xo,r) 



\x-xo\ \x-xo[ 



'The relation l|2.22|l is usually called exact atomic balance. 
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for some < a < 1 and some k < \/3/2, as is obviously seen from the proof. 

We have proved that the widely used kinetic balance method allows to avoid pollution in 
the upper part of the gap for smooth potentials, hence for instance for V — — p * |a;|~^ where 
p > is the distribution of charge for smeared nuclei. However, the kinetic balance method 
does not avoid spectral pollution in the case of pointwise nuclei (Coulomb potential). 

Proof. 



Case (i). We assume that V e ^^(R^) n L°°{m.^) satisfies ([3J8)) . Clearly we have 
sup^ fj-i{ip) < — 1 + sup(V') =: d and m2 = infi^ ^^2{f) > 1 + inf (V^). Hence we necessarily 
have m'2 > > d as requested by Theorem 13.21 Also since V is bounded by assumption, 
(D° + V^)|PL2(K3 £4) = 1 + is bounded, hence p.Sp holds for b — 1. We deduce that 

SpniD'> + V,r,LKB)n{c,l) C [m'2,1). 

Now we claim that > 1. Indeed, let us argue by contradiction and assume that there 
exists a sequence ipn & C^{M.^,C^) such that tpn — ^ in L^, ||(/3„|| — 1 and fi2{'fn) ^ A S 
(c, 1). The number /i2 (</?„) is characterized by the equality 

V\^n\'' + - ^ ^ =M2(¥'n)-l- (3.20) 

{H2{'fin) + l-V)\a-Wipn\^ 



1 

JR3 



Since V is bounded and \\ipn\\ ~ 1 we get 



< (1 - A + 0(1) + ||y||^)(l - A + o(l) + \V\^) / \o ■ V^. 



|2 



2^ 



which proves that {^n) is bounded in 7J^(R^,C^). We deduce that iy9„ ^ in LP(R^,C 
weakly for all 2 < p < 6 and strongly in L^^^(R^, C^) for all 2 < p < 6. Under our assumption 
on V , this shows that lim„^oo / ^I'i^'nP = 0. For n large enough, we thus have 

2 



(P2(93„) + 1- V^)|cr- V(p„P 



^<0 (3.21) 



which is a contradiction since by assumption [i2{!^n) = A + o(l) > d~>V —\. Hence we have 
proved that Spu(i:)° + "P, L^s) n (d, 1) = 0. 

Now we assume sup(V^) > (otherwise there is nothing else to prove since d — — 1) 
and prove that (—1, —1 + sup(V')] C Spu(I?° + V,V, Lkb)- Let xq be a Lebesgue point of 
V, with V{xo) > (hence V{x) > on a neighborhood of xq). Consider a smooth radial 
nonnegative function C which is equal to 1 on the annulus {2 < \x\ < 3} and outside the 
annulus {1 < |a;| < 4}. We define for some fixed 5 > 

(„"^cw.-.o))+^c(^))(; 

where we have chosen the scaling in such a way that the above two functions have a disjoint 
support. We note that 

J \ip^f =5N + 0(n-2), J \a ■ Vip„f = D + 0{Sn-^) 
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where we have introduced N :— J and D — J |VCp. Similarly, we have, using (|2.35|) and 
that V ^ at infinity, 

((1 + V)ipn, ^n) = 6iN + o(l)) + 0{n-^), 

((-1 + V)LKBiPn,LKBiPn) = (-1 + V{xo))D + ©(n'^). 

Hence the matrix of D'^ + in the basis {{(fin, 0), (0, LKBfn)} converges as n — > oo towards 
the following 2x2 matrix: 

( 1 ] 

Eventually we note that ipn WfnW'^ and a ■ V^Pn ic ■ V(p„||~^ ^ 0. Hence, varying 5 and 
ccq, we see that Mi = — 1 + sup(y) and mi < —1 where mi and Mi were defined in p.ip . 
This ends the proof of (i), hy Theorem 13. II 

Case (ii). We will use again Theorem l3.1l More precisely we will show that m2 = — oo < — 1 
and M2 > 1, where m2 and M2 have been defined in (|3.2|) . This time we define 



M^) = (n'^'Cinx)) + (SnY^HiSnx))) (^J^ (3.22) 

where (5 > 4 is a fixed constant (note the above two functions then have a disjoint support). 
Similarly as before, we compute 

{{1 + V)(pn,fn) ^ -n Ci, 

((-1 + V)LKBVn,LKB'fn) = -2D - k(1 + 5)nC2, 

where N and D are defined as above and 

d = / — j— j — dx, 02=/ i— i ax. 

Jr^ m Jr3 \x\ 



.... _ I i~«^*T+r^iv "-((T+r^j 

/ . 1/2 

Let us now choose 5 > 4 large enough such that k^(1 + (5"^)(1 + 5)CiC2 - 2D^ > 0. Then 

det(A„(5)) = (1 + ,5-2)jV£) n ^0{n) (3.23) 

hence det(A„((5)) ^ +00 as n ^ 00. Note that the first eigenvalue fii{(pn) of An{S) satisfies 

C2 

Mi(<^«) < -1 - k(1 + S)n— 
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hence /ii((p„) — cxd as n ^ oo. Therefore we must have ^2{^n) < for n large enough. 
More precisely 

1/2 

"20 '" \{1 + S-^)N^ 
therefore, multiplying by /i2(Vn) a-nd using (|3.23p we deduce that 

m{^n) < TT^n + 0{1) 



, . , C2 f 2D 

p^iWn) > -1 - k(1 + o)n— - n 



which eventually proves that fi2ifn) —» —00. As it is clear that (p„ Iv'ni ^ and a ■ 
V^pn \iy ■ V(/5n|| ^ 0, we have shown that TO2 = —00. 
The proof that M2 > 1 is simpler, it suffices to use 



^„(x)=n-3/2cQ 



whose associated matrix of A reads 



1 - 

'^Nn y N n 



Therefore the result follows from Theorem 13. II □ 



3.2.2. Atomic Balance 

We have proved in the previous section that the kinetic balance method allows to avoid 
spectral pollution in the case of a smooth potential, but that it does not solve the pollution 
issue for a Coulomb potential. In this section we consider another method called atomic 
balance. It consists in taking 



Lab = ■ i-^^) 



(3.24) 



where we recall that we have assumed 2 > sup(l^). Provided that V is smooth enough, we 
can define Lab on the domain D(Lab) = C^{M.'^ \ {0}, C^), in which case Lab satisfies all 
the assumptions of Section [3) Our main result is the following 

Theorem 3.4 (Atomic Balance). Let V be such thatsup{V) < 2, {2~VY'^W e L°°(R^) 
and 

for some < k < VS/2. We also assume that the positive part max(V,0) is in if(IR^) for 
some p > 3 and that lim|j;|^oo Vi^) — 0. Then we have 



Spu{D" + V,r,LAB) = [-1,-1 + supT/]. 

Remark 3.6. We define the operator Lab on D{Lab) = C^{^^ \ {0},C2). Note that 
under our assumptions on V we have that LabD{Lab) is dense in H^(M? ^C"^) for the 
associated Sobolev norm, hence Lab satisfies the properties required in Section \3.1.S\ . The 
above conditions on V are probably far from being optimal. 
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Remark 3.7. The choice of '2' in the definition of Lab is somewhat arbitrary. It can be 
seen that our result still holds true if sup (V) < 1 and Lab is replaced by (6 — V)~^a ■ p for 
some fixed > I. The proof is the same when 9 > 2 but it is slightly more technical when 

i<e <2. 

As we will explain in the proof, a very important tool is the Hardy-type inequality: 

This inequality was obtained in [12' by using a min-max characterization of the first eigen- 
value of ~ica- V+c^P-iy/lxl Indeed ((3:251) is an equality when tp is equal to the upper spinor 
of the eigenfunction corresponding to the first eigenvalue in (—1, 1) of —ica ■W + c'^/3—i'/\x\. 
The inequality (|3.25p was then proved by a direct analytical method in [TT|. Introducing 
m = c(l -|- -y/l — (ly/c)-^) and k = v j c we can rewrite (I3.25|) in the following form 



^' dx + rt-i / \ip[x)\^dx>K / dx. (3.26) 



We now provide the proof of Theorem 



Proof. Let us first prove that when sup(l/) > 0, then we have (—1,-1 + supV^] C 
Spu(£'° -I- V.,V tLab)- The proof is indeed the same as that of Theorem 13.31 we define for 
some fixed 5 > 

I \ ( in/-t I w , ^^^'^ ^fx-xo\\ (\ 

ipn\x)=\n ^ C \n{x - Xq)) + , C 



(4n)3/2^ V 4n 77 VO 

where Xq is a Lebesgue point of V such that < V{xq) < 2. One can prove that the matrix 
of + y in {{fn, 0), (0, Lab'Pu)} converges as n ^ oo towards the following 2x2 matrix: 

i^r ] 

-l + V{xo)l 




Hence we have again, by Theorem 13. 1[ (— 1, — 1 + sup V] C Spu(_D° + ¥,1^, Lab)- 

The second part consists in proving that there is no spectral pollution above — l-l-sup(y). 
As a first illustration of the usefulness of the Hardy- type inequality p.26p . we start by 
proving the following 

Lemma 3.2. We have 



m" = inf n2(ip) > 1 - 2- (3.27) 

Remark 3.8. We note that the right hand side of (j3.27p is always > 1/3 when < n < 
VS/2, and it converges to 1 as k 0, as it should be. 

Proof. The number iJ,2{f) is the largest solution of the equation 



/ r |CT-Vy(a)|^ , 



2-V{x} 



2 



il + Vix)mxr+ ^ a+^-y(.))|..v^(.)P , =m/ l^ix^dx. (3.28) 

Jr3 (2-V(x))^ "-^ 
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Clearly we must always have 

Jr3 (2-y(a;))2"^ 

Let be ndf) < ^ < 1. We estimate: 



2 

2, 



> (1-2^^^^-^) / (3.29) 

y 1 + vi - y jR3 

where in the last line we have used ()3.26p and the fact that K|a:;|^^ + V{x) > 0. From this 
we deduce that 



2 



/ 1 — vT — K 

fi2[(p) > max I 1 - 2— — ^ , Mc(<(5) 



This ends the proof of Lemma 13.21 □ 
The next step is to prove that property (|3.8p is satisfied. 



Lemma 3.3. Property (|3.8p holds true for b = oo: if {(pn} C C5"(R'^,C'^) is such that 
(fin in and ^2{fn) ^ i < oo, then J^^ l/|iy9„p — > 0. 

Proof. Note that necessarily i > 1/3 by Lemma [3.21 hence £ must be finite. We use the 
estimate p.29p . with ^ = H2{Vn) to get 

0> / il + Vix)-ii2{'Pnmnix)\^dx+ [ t^l^f^l^dx. (3.30) 



2^R 



Now we write 



>(1-^)/ + ^d.-/ Mx)\^dx 

7r3 ^ + R "'R^ FI JM3 

(3.31) 

where we have used (j3.26p with m <-> 2/k and k 1. We deduce that 

+ Vix)]\ipUx)\'dx+{l-K^) f \^-l^£ll^dx<fl2{iPn) f \iPn\'. (3.32) 



/r3 V fI / 2 r 

Using that ^2(Vn) ^ that V > — and ^ we deduce that 

/• |a-V^„(x)|% 
um / ax = 0. 

"-^"^ Jr3 2 + R 

Using again (|3.3ip with ip ~ ipn we finally get the result. □ 
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We will now prove the following 
Lemma 3.4. We have m'2 > 1 where m'2 was defined in (|3.7I 



Proof. Consider a sequence {^Pn} C C^(M'^,C'^) such that ||<y9„|| — 1 and (y9„ 0. We 
will argue by contradiction and suppose that, up to a subsequence, p.2(}Pn) ^ £ & [1/3, !)• 
Similarly as in the proof of Lemma |3.3[ {(/?„} must satisfy p.32|) . from which we infer that 



2+ A 

\x\ 

hence {ipn} is bounded in H^. Therefore, up to a subsequence we may assume that ipn 
strongly in Lf^ ^(M^) for 2 < p < 6. 

Let us now fix a smooth partition of unity Co + + ~ 1 where each is > 0, = 1 
on the ball B{0,r) and = outside the ball 5(0, 2r), 6 = 1 outside the ball 5(0, 2i?) 
and ^2 = in the ball B{0, R). We fix R large enough such that 

V|x| > R, \V{x)\ < ^ 

and r small enough such that 

e 

m — e < — 
- 2r 

where e is a fixed constant chosen such that 1 — I ~ e/3 > (1 — £)/3 and k + e < \/3/2. 
Next we use the (pointwise) IMS formula 



|V^(x)|2 = ^ |V(6^)(2:)|2 - \^{x)\' ^ |Ve.(x)|^ 

i=0 1=0 

and p.30p to infer, denoting (pl^ ipn^i and 77 = I]Lo l^6(2;)P, 

jzlf (i + nx)-/i2(^„))K(x)pdx+ / ^i^^^dx") 



< / (3.33) 



Next we note that for n large enough, by our definition of R, 

(1 + Vix) - f,2iiPn))\^lix)\^dx > i-^ W^lf . (3.34) 
1 o 

Similarly we have by definition of r and e (using that yj^j has its support in the ball -6(0, 2r)) 



where for the last inequality we have used (|3.26p and k + e < \/3/2. Using again that 
V > —k\x\~^, we infer the lower bound, for n large enough, 

/ {l + V{x)-^,2{'Pn)M{x)\'dx+ f t:l^pldx>^y'if. (3.35) 
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Inserting ()3.34|) and ()3.35|) in ()3.33|) . we obtain 

l^-^/ll 2||2 , II 0||2\^ f V{x)\'fn{x}\'^ 



f\\ 2||2 , II 0||2\^ / V[X ^^(X \ , \\rr^ II II -t ||2 

[\\^n\\ +hn\\ )< 2+^ + \\Vlr<\x\<2R\\L^\\^ntr<\x\<2R\\L^ ■ 

Using the strong local convergence of (p„, we finally deduce that lini„^oo H^nll — 
lim„_,tx) — which is a contradiction with = 1. □ 



The conclusion follows from Theorem 13.21 (ii). This ends the proof of Theorem 13.41 □ 
3.2.3. Dual Kinetic Balance 

Let us now study the method which was introduced in I32j , based this time on the splitting 
of the Hilbert space induced by the projector Vc defined in (|2.29p . We have seen in Theorem 
12. 51 that pollution might occur when e is not small enough. We prove below that introducing 
a balance as proposed in [32j does not in general decrease the polluted spectrum. 
Let us introduce the following operator 



J 



defined on VL'^{^^ X"^) with values in {I - V)L'^{M? .C"^). Next we introduce the following 
balance operator [32] 



Ldkb^U.JU, (3.36) 



which is an isometry defined on V^L'^iM.^ ,C'^) with values in {l~re)L'^{M.^ , C^)- A calculation 
shows that, like in [32^, formulas (24) and (25), 



L 



DKB 



As before we may define Ldkb on C = U,C^{M?, C^). 

Theorem 3.5 (Dual Kinetic Balance). Assume that the real function V satisfies the 
same assumptions as in Theorem \2.4\ Assume also that Ve and Ldkb cLfe defined as in 
(|2.29p and p.36p for some < e < 1. Then one has 



Spu(DO + y, P„ Ldkb) = Spu(i?" + V,V,) 

r 2 11 r r 2 

— 1 , min < h 1 + sup , 1 > U max i ~ 1 j 1 + ^ 

Proof. We will use Theorem 13.11 Consider a radial function C e C^(R'^,M) and intro- 
duce the following functions: ipi := (C,0) and Lp'^ := (ct ■ p)/\p\ipi £ ns>o-ff'*(R^, C^)- 
We define similarly as in the proof of Theorem 12.51 (pn{x) = n^^^ipi{n(x — xq)) and 
(p'nix) = n^^^(f[{n{x — xa)), where xq is a fixed Lebesgue point of V. We note that 
(p'„ :— (a ■p)/\p\ipn. Also, using that C is radial, we get for any real function /: 

(/(|p|)^„, p'J = (fin\p\)ip^,ip[) = I Lo.dLo r \a\p\)\^f{n\p\)\p\''d\p\ = 0. (3.37) 

Js^ Jo 
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A simple calculation shows that the 2x2 matrix of D'-* + V in the basis 
([/e(v5„,0) , LoKBUei^Pn-.^)) reads 



M„ = 



where 



All = 1 + -^=^==V ^ 



ecr ■ p / 2 



2 V e 



2 + y 



6(7 • p 



= -1 



1 -.V- ' 



ecr • P / 2 



A12 = (A21)* = -^—^^(a-p)+e-^===[V,a-p]- 



ecr • p 

1 



We infer from p.37p that 



for every n. Also we have 



It is then easy to see that 



2e-l + e^\p\\ 



lim 



= 0. 



hm M„ = ^ n 



+ i + y(xo) 

Note that i£)A's(<^ni 0) — ^ since Ldj^b is an isometry. The result follows from Theorem 
13. H by varying xq. □ 
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